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Abstract 

We have studied the U(l) gauge field theory with Villain (peri- 
odic Gaussian) action on spherelike lattices. The effective size of the 
^f) . systems studied ranges from 6 to 16. We do not observe any 2-state 

I signal in the distribution function of the plaquette expectation value 

' at the deconfining phase transition. The observed finite-size scaling 

\^ • behavior is consistent with a second order phase transition. The ob- 

i tained value of the critical exponent is u = 0.366(12) and thus neither 

Gaussian (i^ = 0.5) nor discontinuous {u = 0.25) type, indicating a 
nontrivial continuum limit. 
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1 Introduction 



In this work we study the Villain formulation of the U(l) pure gauge 
theory on 4D lattices with spherelike topology. We were motivated by recent 
results obtained in the theory with Wilson action (with and without a charge 
two coupling 7) on spherelike lattices |21 El El- It turned out that on such 
lattices there are no 2-state signals even on the Wilson line (7 = 0) [3]. 
For hypercubic lattices with periodic boundary conditions (the usual torus 
geometry) one finds such signals El E] even at sizable negative values 
of 7 0; under certain assumption a tricritical point (TCP) was predicted 
at7 = — 0.11|ni. It has been argued |2j that the disturbing effects may be 
due to the interplay of periodic boundary condition with the important role 
of the topological monopole excitations in the transition. This issue is still 
unsettled, however. In case the two-state signal persists in an infinite- volume 
limit, the transition is first order; if the signal is spurious, it prevents a careful 
finite-size scaling (FSS) analysis of the critical transition. Up to date taking 
a continuum limit has not been possible due to this problem. 

We hope that in this situation simulations with the Villain action help to 
clarify some of the problems mentioned above. This formulation is of certain 
interest because it lends itself to theoretical analysis. In this form of the the- 
ory the partition function can be decomposed into Gaussian fluctuations and 
monopole excitations |n|. There are further relationships to the noncompact 
U(l) Higgs model in the limit of large negative squared bare mass (frozen 4D 
superconductor) [TUl E] and to an effective string theory equivalent to that 
model OCn]. The action obeys reflection positivity. 

The leading terms of a character expansion of the Villain action shows 
that the phase transition (PT) is near Pw = 1-16 (we denote by Pw the 
coupling in the Wilson action) and 7 = —0.22. Although this value of 7 is 
below the value for a conjectured TCP |^ metastability signals were observed 
(for torus b.c.) here, too ^lE]. The PT lies in a hypersurface of PTs for 
the other actions mentioned. It is therefore of interest, whether the two-state 
signals vanish for spherelike lattice for the Villain action as well, and, if this 
is the case, whether the critical behavior is in the same universality class. 

We proceed in the following way: In section 2 we discuss the details and 
implementation of the Villain action, in section 3 the lattices with spherelike 
topology are introduced and the general strategy of the simulation and the 
FSS analysis are presented. Section 4 summarizes our conclusions. 
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2 The Villain action 

The Villain or heat-kernel action is defined through the Boltzmann factor 
per plaquette 

exp(-5p(Af/p)) = ir(l/Af/p), (1) 
and is the solution of the diffusion equation in group space 

AK{t,U)^-f^Kit,U), (2) 

with the boundary conditions 

K{0, U) = 5(1, [/), K{oo, U) = 1. (3) 

It connects the constant distribution at strong coupling {t = oo) with a 
distribution peaked at the group unit element in the weak coupling limit 
(t^O). 

The heat equation is implemented by introducing the metric tensor and 
the Laplace-Beltrami operator on the group manifold. In the case of the U(l) 
gauge group the solution can be written as 



exp(-5'p) = 



(4) 



where ^p is the angle of the plaquette variable Up — cos9p and given as a 
sum over all link angles O^^n G [— tt, tt] 

In the simulation Sp has been calculated by interpolating pretabulated val- 
ues. 

The total action is defined through 

exp (Sm = IJexp (-Sp((3, Up)) (6) 
p 

and the observables (quantum averages) are given as expectation values with 
regard to the Gibbs measure Y[x,nd'Ux,f^exp {—S{P)). Useful observables are 
the internal energy 
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and the specific-heat 



This form differs from the usual definition due to the non-hnear dependence 
of Sp on p. Note, that /? is different from the couphng used in the 
Wilson-action. 

In our calculations we measure the expectation value of the Wilson pla- 
quette variable 

E={Y,wpReUp) (9) 
p 

and its distributions (the weight factors are due to correction factor for the 
spherical lattice shape and are discussed below). For the investigation of the 
PT both sets of observables are equally suited; we find the later set more 
useful, however. The histogram analysis along the lines of Ferrenberg and 
Swendsen ^H] is implemented in that variable E and the analytic continua- 
tion therefore is in the conjugate coupling variable f3w From the histograms 
in E corresponding higher order cumulants are determined, 

cv{(3,L) = ^{{E-{E)r), (10) 

Vclb[P,L) - ^^^p , (11) 

U(3L) - (12) 

(For simplicity we call cy specific-heat, too.) The positions and values of 
their respective extrema are used for the FSS analysis. 



3 Numerical simulation and results 

We have performed simulation on the following lattices with spherelike topol- 
ogy: 

• SH[N], the surface of 5D hypercubes of size N^; 

• S[N], with the topology of SH[N] but introducing weight factors for the 
plaquette action: S = J^p'i^pSp, correcting for spherical shape, i.e. 
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distributing the curvature over the lattice (details about the geometry 
and these factors are given in Q). 

Most of our simulations were done on S[N] lattices with N =4, 6, 8, 
10; these lattices have roughly the same number of variables as 6^, 9^, 12^, 
16^ hypercubic lattices. In average we have performed about half a million 
sweeps for each lattice size, using a 3-hit Metropolis algorithm. We have 
measured the values and distribution functions (histograms) for Ehk and E 
and found no indication of double peaks. Unfortunately a multihistogram 
analysis combining various histograms determined at different values of (3 is 
not applicable here because of the specific form of the action: Sp does not 
depend linearly on f5. However, we are able to do the analysis based on single 
histograms in E in the vicinity of the PT and extrapolate into the "Wilson 
direction" (3w based on the individual histogram, 

ifiE)) = 4 E h{E, P)f{E)e-^^\ A^Y. h{E, /3)e-^-^, (13) 

^ E E 

where h{E,f3) denotes the histogram entry at E. The number of bins was 
taken sufficiently large (typically > 2000) to exclude systematic errors. 

For each lattice size we try to simulate as close as possible to the peak 
position of the specific-heat. For this aim we first determined the cumulants 
as functions of jSw from histograms at various values of (3. This preparatory 
analysis was done at the statistics of lOOK sweeps. From this we infer a 
relationship between Py[r{peak) and /5 and determine the presumed value of 
the corresponding pseudocritical coupling /3o. These are given in Table 1. 
Then we performed long runs (about 500K sweeps) at these couplings. The 
final results are histograms determined almost on top of the specific-heat 
peak. The individual histogram analysis allows us to extrapolate away from 
this point in coupling space into the direction and the table also shows 
this distance to the actual peak position of the specific- heat. This provides 
further (small) corrections. 

For our analysis we use the values of the cumulants at their (extrapolated) 
maxima or minima, which agree within the errors with the values measured 
directly at Pq. 

The final histograms are similar to those determined for the Wilson action 
|3] and show no 2-state signal. (Note that these are individual histograms and 
not the result of a combination.) We proceed to study the FSS of the extrema 
value of the specific-heat (maximum) and Binder cumulants Vclb and U4 
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S[N] 




A/3h' 


S[4] 


0.63687 


-0.00029(46) 


S[6] 


0.64550 


0.00060(28) 


S[8] 


0.64766 


0.00067(10) 


S[10] 


0.64867 


0.00017(14) 



Table 1: Couplings (3o where the long runs were performed; we also give the 
extrapolation distance to the peak position of the specific-heat in direction 
of Pw (as determined a posteriori from the histogram analysis). 

(minima, cf. [21 E] for niore details on the definitions and FSS properties) 
to check for consistent critical (i.e. second order ) FSS behavior, which in 
leading order may be parameterized 

cv,maAL) ^ a, + KL^/\ (14) 
VcLB,n.^niL) ^ (a, + 6,L"/'^)L-^ (15) 

U^,rmniL) ~ a„ + 6,^""/^ (16) 

In our analysis we assume validity of the hyperscaling relation a = 2 — Du 
[D = 4 in our case). The length scale L = V*, where the volume V = 
I Z^p'W^P) roughly proportional to the number of sites (in fact, this definition 
gives the number of sites in the situation of the usual hypercubic lattice with 
torus boundary conditions). 

A log-log plot of the peak value of the specific-heat as a function of L 
is shown in Fig. la. From the fit we obtain the value ajv = 1.479(135) 
and V = 0.365(9) {x^/d.f. = 1.3), which is far from what one expects in the 
case of a first order PT, where a/u = A and u = 1/D = 0.25. This value is 
also compatible with the values obtained in 00]. The cumulant VcLB,min 
is shown in Fig lb; its FSS behavior is in good agreement with that of the 
specific-heat. The fit due to ^ leads to u = 0.359(10) (x^/d.f. = 0.8). 
The cumulant Ui^min is compatible with the expected FSS; it approaches 
a constant ~ 2.75(3) but has too large errors to determine the nonleading 
scaling term, although it is consistent with the values of v determined from 
the other quantities. 

Let us discuss briefly the FSS of the pseudocritical couplings (as sum- 
marized in Table 1). Generally the FSS of pseudocritical couplings may be 
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Figure 1: The logarithm of the extrema values of (a) the specific- heat and 
(b) the Challa-Landau-Binder Vclb cumulant vs InL; the fits represent the 
leading FSS behavior given in ()14|1 . 

written in the form 

PiV) = Pioo) + cL-\ (17) 

where A is the so-called shift-exponent. For many models X = l/u, but this is 
not necessarily so in general (see the discussion in a recent study of the 
2D Ising model on lattices with spherelike topology yielded X ^ l/u [H^). 
Therefore the FSS behavior of the pseudocritical coupling is not necessarily 
suitable for an independent determination of u. A fit to the values in Table 
1 gives A = 2.71(26) - consistent with X = l/u within the large error bars 
- and P{oo) = 0.6496(3). Note, that due to the geometry correction factors 
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this latter value does not have to agree with the corresponding value for torus 
topology. A more detailed discussion in a similar context can be found in 
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Figure 2: The imaginary part of the closest Fisher zero as the function of 
the size L. 

Another efficient way to determine the exponent u lies in the FSS of the 
imaginary part of the (Fisher) zeroes zq of the partition function (in the 
complex coupling plane) closest to the real axis. The corresponding plot is 
shown on Fig.2; the expected scaling is Im2;o oc L^^^'^. From the fit one 
obtains z/ = 0.366(12), using all lattice sizes {x^/d.f. = 0.5). This is in 
excellent agreement with the values obtained from the cumulants. 

We have also performed simulations on SH[N] lattices with L = 4, 6, 8. 
Since here the statistics was only 20% of that for the S[N] lattices we do not 
quote errors and did not analyze these data in more detail. The main result 
is that also for SH[N] we do not see any 2-state signal. The pseudocritical 
couplings Po are 0.64, 0.64375, 0.645 for corresponding lattice sizes. The 
peak values of the specific-heat are 1.48, 2.77, 3.32 correspondingly. All 
these results are close to those obtained on S[N] lattices. 




8 



4 Conclusion 



We have studied the compact U(l) gauge theory with the Villain action on 
4D lattices with spherelike topology. We do not observe any metastability 
(i.e. 2-state signal) in the distribution function of the internal energy. The 
FSS for various quantities is consistent with critical behavior. The obtained 
value of the critical exponent v is compatible with neither the Gaussian value 
0.5 nor the "discontinuity" value 0.25. In our determination for the Villain 
action we get values in the range 0.36-0.37, a smaller range of values, but in 
perfect agreement with those obtained in ^ with mixed action on spherical 
lattices. Note that these actions are at quite separate points in the space of 
plaquette actions given by a character expansion. 

These results indicate that below the Wilson line the deconfining PT is 
of second order and both the mixed action and the Villain action belong to 
one universality class of models with a nontrivial continuum limit. 
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